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I. INTRODUCTION 

The motion of self gravitating systems and the associated gravitational waves have been 
among the major research interest in general relativity. Furthermore, understanding systems 
with strong internal gravity such as neutron star binaries are important also in astrophysics, 
since they are thought as a possible candidate for gamma ray bursts [I[] as well as the most 
likely sources of gravitational waves. Laser interferometric detectors are under construction 
for detecting gravitational waves (LIGO §, VIRGO §, TAMA |[ and GEO600 §). It is 
expected that these will become new eyes looking into the universe and bring great impacts 
on astrophysics. For this to become realistic, reliable comparisons between theoretical pre- 
dictions and the observations are definitely necessary. Namely we have to prepare sufficiently 
accurate theoretical templates for waveforms for likely sources including inspiralling neutron 
star binaries. This is possible only if we know their precise equation of motion 1^-^. Since 
the Einstein equation is highly nonlinear system as is well-known, the use of some kind of 
approximation is unavoidable to describe a dynamical system such as a neutron stars binary. 

Among various approximation schemes, the post-Newtonian approximation is most pop- 



ular for treating nonlinear gravitational interactions |10], and the equation of motion has 



been derived successfully up to the 2.5PN order |Tl|-p3|. At higher orders, however, an ex- 
pansion in power series seems to break down (e. g. [T^). Hence, it is very interesting to see 
what form the equation of motion takes at higher orders. In recent, 3PN, 3.5PN and higher 



orders are being tackled [0,^. Furthermore, some methods have been developed to treat 



the matter of binaries and internal gravity of them. For instance, in the post-Minkowskian 



approximation [|Tl|)|ll] j the self gravity causes divergent behavior at the position of a point 



source, so that regularization schemes have been used [|T^,|r2]. Other method assumed that 



the neutron star is modeled by the spherical distribution of the perfect fluid pO| , p~^ . In these 
derivations where the standard post-Newtonian approximation is used, however, the weak 
field is assumed everywhere or everywhere except at the singular point. As a result, strictly 
speaking, the obtained equation of motion is applicable only to a object with weak internal 



gravity or to a singular source whose relation to the real system is not clear. It is then 
necessary to extrapolate the equivalence principle to the extended objects with strong self 
gravity, in order to apply the post-Newtonian equation of motion derived from the standard 
approach to neutron star binaries. Since such an extrapolation is not trivial at all, it must 
be shown within the framework of general relativity if it works. 

The equation of motion for bodies with strong self gravity has been studied by various 
authors; D'Eath developed an approach based on the asymptotic matching |^. However, 
this approach is based on a perturbation of black holes. Hence, there remain several problems 
to be solved in this scheme; how can we take into account the difference between neutron 
stars and black holes? How can we go to nonlinear perturbations? (Recall that the nonlinear 
perturbation is required for equal- mass binaries even at the 1 PN order.) Thorne and Hartle's 



derivation is also using the linear order perturbation ||22|. However, there seem a great of 
difficulties in pushing their method forward to obtain a higher order equation of motion. 

In the present paper, we take other approach, the so-called point particle limit in order to 
treat the strong self gravity where the size of the component star is reduced to a point keeping 
the strength of the internal gravity pB|,|2^. The method has been used to derive only the 



equation of motion at the Newtonian order and the quadrupole radiation reaction at 2.5PN 
order for compact binaries. Here and henceforth, by the post-Newtonian approximation, 
we mean only the smallness of the ratio of the orbital velocity to the light velocity and the 
weakness of the gravitational field outside the compact objects, but not the weakness of the 
gravitational field inside these. As discussed above, we need to know higher post-Newtonian 
orders up to, say, 4th post-Newtonian order to have an accurate prediction of the emitted 
waveform p5[. This is what we are aiming at. Namely, we wish to derive the equation of 



motion applicable to systems with strong internal gravity such as compact binaries up to a 
sufficiently higher post-Newtonian orders using the point particle limit. As the first step, 
here we derive the first and (partly) second post-Newtonian equation of motion for compact 
binaries taking fully strong internal gravity into account. For the sake of the completeness, 
we shall show some details of the approach which has been proposed by one of the present 



author p4|. Our main contributions in this paper are to derive the first post-Newtonian 
equation of motion, the lowest spin-orbit couphng force, the lowest quadrupole-orbit coupling 
force for relativistic compact binaries and a general momentum velocity relation which was 



not discussed in [Q. 

This paper is organized as follows. Section II presents a general framework for the 
derivation; the equation of motion expressed with surface integrals and the point particle 
limit. In section III, we review the Newtonian-like equation of motion for compact binaries 
to show our basic calculation. We derive the IPN equation of motion in section IV, namely 
with some remarks on the definition of the mass and momentum in this paper and take 
account of effects of the quadrupole moment and the spin of compact objects in section V. 
Section VI is devoted to the conclusion. In appendix A we will give some comments on the 
derivation of the components of the metric. Some formulae are given in the appendix B. 

Greek indices take from to 3, and Latin indices from 1 to 3. 

II. GENERAL FORMALISM 

A. Newtonian scalings and point particle limit 

It is expected that the orbital motion of neutron star binaries in the beginning of the 
inspiralling phase are governed mainly by Newtonian gravity so that it is natural to formulate 
our approach based on the post-Newtonian approximation. For our purpose it turns out 



convenient to introduce a Newtonian dynamical time as |16 



Then, the coordinate velocity is 



T = et. (2.1) 



^ = e^, (2.2) 

at dr 



so that the velocity scales as e. Note that events with different e but at the same Newtonian 
dynamical time r are at (roughly) the same phase in the orbit. This motivates us to define 



the post-Newtonian sequence by letting e go to zero with r staying constant. We call this 
limit, e — > 0, the Newtonian limit. In the binary whose mass and separation are denoted by 
M and L respectively, the balance of Newtonian forces is written as 



|2 



dt 



r^u 



M/L. (2.3) 



Hence, we find that M scales as e^ when we fix the separation L. 

In order to treat strong internal gravity we introduce the point particle limit in the 
post-Newtonian sequence. Namely, denoting a size of a neutron star by R, we expect M/R 
is of order of unity. To be so, R must scale as e^, which means that the density is of 
order e~^. Hence, we define the body zones Ba and the body zone coordinates as Ba = 
{x'^Wx'' — -2^(t)| < c-Ra} for some Ra and 

a^ = e-\x'-z'A{T)), (2.4) 

so that in the body zone coordinate, a size of a neutron star remains constant even as e ^ 0. 
Here the underlined indices refer to the body zone coordinate and Za{t), A = 1,2, are the 
near zone coordinates of the representative points of two bodies. There are three purposes 
to define the body zone and the body zone coordinate. First, it is useful to define physical 
quantities such as mass, spin with the body zone coordinate, since in this coordinate the 
body is nearly isolated. Second, as can be seen later, we define a gravitational force on 
the object A as a momentum flow per unit time through the surface of the body zone of 
the object A. Thus we only need to know the metric at the boundary of the body zone 
where the assumption of weak gravity is valid. To deal with strong internal gravity, this is 
a great advantage since we do not have to care about the internal strong gravity. Third, the 
multipole moments of the body are reduced to the higher order terms relative to monopole 
terms thanks to the scaling of the body zone coordinate, as will be seen later. This is a nice 
feature because in the inspiralling phase it is a good approximation that compact objects 
are expressed as point particles with their multipoles as higher corrections. 

In this paper we adopt the initial value approach as in PB|,p3 and we choose two sta- 



tionary uniformally rotating fluids with spatially compact support as the initial data for 



the matter. Also, we assume that both objects rotate slowly, that is, the velocities of the 
spinning motion of the bodies scale as e. It is easy to incorporate rapidly rotating stars into 
our formalism. The result is that a coefficient e^ appears instead of the coefficient e^ in the 
equation ( |5.14| ). We note that, with this scaling, the spin-orbit coupling force appears at 



the 2PN order, which is different from the usual result [|r^,^,^. The initial data of the 
stress-energy tensor T^'^ of the body A is supposed to scale like 

T7 = e^T^* ~ e"^ (2.5) 

tJ = e-^T'2 ~ e-^ (2.6) 

r^=e-^r/ ^t-\ (2.7) 

Note that as a slowly rotating fluid, the object is supported by the pressure. Also, note that 
the assumption that the densities of the fluids scale as e~^ are applied to the components in 
the near zone coordinate, not in the body zone coordinate. The scalings of the stress energy 
tensor of bodies on the initial hypersurface can be found in the Table |. 

Before ending up this subsection, we mention briefly the coordinate transformation from 
the near zone coordinate to the Fermi normal coordinate of the extended body which is 
needed to derive the spin-orbital coupling force in the same form as the previous works 
11,11111 • 



These transformation may be obtained order by order, and at the lowest order, we identify 
this transformation as the Galilei transformation, so that we have 

T^^ = T7, (2.8) 

T- = e^TP + TYv\, (2.9) 

Ti^ = Tr^>i + ^e'v^Tt + e'Tl (2.10) 

where the subscript N means that these are components of the stress energy tensor in the 

near zone coordinate. These equations are enough for obtaining the IPN equation of motion. 

It is difficult, however, to obtain the transformation even at the next order when we 

consider a body with strong internal gravity (for the case of weak internal gravity, see e.g. 

6 



29| , |30[] ). Fortunately, it turns out in section V that we do not need to know an explicit 



form of the transformation to derive the lowest spin-orbital coupling force. 

B. A general form of the equation of motion 

As discussed above we need to know the metric at the boundary of the body zone where 
it is assumed that the metric deviates slightly from the flat metric. Thus we define the small 
deviation from the flat metric rj'^'^ as 

h>"' = 7]^"^ - ^/^g^", (2.11) 

where we denote the determinant of the metric g^p by g. 
We choose the harmonic condition on the metric as 

^1^",^ = ^^ (2.12) 

where the comma denotes the partial derivative. Then, the Einstein's equation is rewritten 



as 



UW" = -167^A^^ (2.13) 



where we defined 



a = r]^''d,d,, (2.14) 

A''" = e"" + X'"""'^ ,al3, (2.15) 

e^'^ = (-(?)(r''^ + tfj, (2.16) 

^M^a/3 _ J_^lau-^I3,^ _ l^Pll^^y (2.17) 

167r 
Here, T'^'^ and t^*^ denote the stress-energy tensor and the Landau-Lifshitz pseudotensor 
3T| , respectively. The conservation law is expressed as 



A'^^,^ = 0. (2.18) 

The reduced Einstein's equation is solved formally as 



JciT,x>']e) \x — y\ 

where C(r, x'^;e) means the past hght cone emanating from the event {t,x^). To clarify 
that A^'^ appearing in the Eq. ( p.l9| ) are components in the near zone coordinate (r, x^), we 
added a subscript A^ to A^*^. This equation must be solved iteratively to obtain an explicit 
form of the metric. We have ignored the homogeneous solution, since it is irrelevant to the 



dynamics of the binary at the order considered in this paper p7 . 

Let us define the four momentum of the object A as the integral over the body zone Ba 
hke 

P^{T)=e'f rfVA]v^^,4;e). (2.20) 

JBa 

Note that x^'^"^ does not play any role up to the 2PN order. We can understand this from 
the fact that the lowest order of the fields h'^'^ is e^, which is seen in the next section, and 
the volume integrals of x^'^°'^ ,ai3 can be transformed into the surface integrals. 

It should be noted that P^ depends only on the time coordinate. We call the n = t and 
i component of P4 simply the energy and the momentum of the object A respectively. Note 
that we have not specified z\{t) yet. We define the derivative which does not change the 
region Ba as 

n ;^ p> 

(2.21) 





D 
Dt 


9 k d 

~ 1 Va 

dr dxk ' 


where we defined v\ as 







v'a = 4(r), (2.22) 

and the dot denotes a temporal derivative. Note that the differential operator D/Dt differs 
essentially from the usual Lagrange derivative d/dr. But since the body zones for each body 
remain unchanged (in the near zone coordinate sense), these two operators are effectively 
the same when they act on the integrals over the body zones. 
Here, we define the near zone dipole moment of the object as 

8 



D'Air) 



Iba 



(2.23) 



To define the center of tlie mass z\{t) of the object A it is required that D\ takes a certain 
value. At the Newtonian and post-Newtonian order, we will require that D\ vanishes. But 



to derive the spin-orbit coupling force in the usual form |ll3| , |2^ , |28|| , care must be taken when 
we define the center of the mass as can be seen later. By taking the temporal derivative of 
Dj^, we obtain 

,dD\ 



P\ = Plv\ + Q\ + e' 



dr 



(2.24) 



where we used the derivative defined by Eq. ( |2.21|) . Here, we defined Q\ as 
Q\ = e-^ <f dSkAt{x' - z\{t)} - e-\\ <f dSkA}j{x' - z\{t)], 



(2.25) 



where OBa denotes the surface of the sphere Ba- Equation (p.24|) is thought as the relation 
between the velocity and momentum if we make a certain choice of the near zone dipole 
moment D\, that is, we define the center of mass of the object A. It is worthwhile to 
mention that the velocity and the momentum are not proportional to each other, but Q\ 
appears as its correction to be considered. Actually, Q\ contributes to the equation of 
motion as shown in the section V. 



Using the conservation law, Eq. ( 2.18| ), we obtain 

Joba Joba 



d_ 
dr 



N ' 



(2.26) 



IdBA JOB A 

where we used Eq. ( [2.211 ). This is an evolution equation for the four momentum P^{t). 

In order to obtain the equation of motion, an evolution equation for the momentum must 
be changed into an evolution equation for the velocity. Inserting Eq. ( p.24| ) into the i-th 



component of Eq. (|2.26|) , we obtain the general form of the equation of motion as 



^4"^ 



dBA JdBA 



dSkA]^ 



+e-\\U dSkAt^-v'xi dSkA 

^JoBa JdBA 



d_ 

dr 



-Q\ 



IdBA 

e 



dr^ ■ 
9 



(2.27) 



The above set of our basic equations ( p.25|) and (p.27|) consist only of the surface integrals 
when the representative point of the body A, z\ or equivalently, the value of the near 
zone dipole moment D\ is specified. Therefore, all we must do is to evaluate these surface 
integrals perturbatively. This is a great advantage of our formalism, since our evaluations 
can be done outside the strong self gravitational field. 

III. NEWTONIAN EQUATION OF MOTION 

From the point of view of the point particle limit, the Newtonian order of the equation 
of motion for compact objects has been examined, where Q\ was not considered in P^ . 



However, since Q\ is newly introduced in this paper, we rederive the Newtonian equation of 
motion in this section. We also intend the completeness and to show our basic calculation. 
At the Newtonian order, we use equations (|2.8| ), (|2.9| ) and (|2.1(]|) . The near zone dipole 
moment becomes 

D\{t) =e^ f d'aAAYir, a\- e)a\. (3.1) 

JBa 

So, we define the center of mass of the object A as z\, by requiring that the near zone dipole 
moment vanishes 

D\{t) = 0. (3.2) 



Then metric up to 0(e) becomes ||23[ 



/,- = 4e^ V ^ + 0(e^), (3.3) 

A=l,2 ^A 

/," = 4e4 ^ :^+4e4 ^ !!M!^ + o(e^), (3.4) 

A=l,2 ^A A=l,2 ^A 

A=l,2 ^A ^=1^2 ^^ A=l,2 ^^ 



-8e^ 



17117712 



, .4 V- "'A i , a OO^K'^2 i J 

A=l,2 ^A ^12'^ 



8\^h - \^''hl 



{ni2 - niy'inu + ^2)'' 



52 
+ 0(6^), (3.5) 

10 



where r^ = x* — z\, ta = |^a|, r\2 = z\ — Z2 and S = ri + r2 + r^- We defined 



m^ = lime^ / d^aA-^V{T,a-A',e), (3.6) 

^^0 Jba 



IBa 



J^ = linie / d aAA'2{T,a-2;e), (3.7) 

Z'i=\ime'' [ rfVAf(r,a|;e). (3.8) 

The r/iA is the ADM mass that the body A would have if the body A were isolated. By 
choosing the center of mass, J^ vanishes. Furthermore, we assume that the fluids of the 
two bodies are (quasi-) stationary so that Z^ are of higher order, e^ [^. Note that at the 



Newtonian order, we have PJ = tua by virtue of Eq. (|2.8|) . We also note that the temporal 
derivative of the PJ vanishes at the Newtonian order from Eqs. (pl|) and (|B^ ). To obtain 
the Newtonian equation of motion, we need only the field h'^'^ of order e^. 
Now, to evaluate the surface integrals we take 

fi = eRifii, (3.9) 

^2 = ri2 + eRiHi, (3.10) 

^2=^1-^2, (3.11) 

where rTi denotes the spatial unit vector emanating from zi and eRA is the radius of the 
sphere Ba- The pseudotensor at this order is 

i4i-9)tfL] = ^{i^)h^^'\^)h^^'' - ls'\^)h^^'i^)h^\i)^ (3-12) 

where (n) means the order of e". We evaluate 



£_ <<5.„|(-9)«?,i = ^^{s'A - 1^'%.) E^ E^i, "sj 



mATTLBr^r^ 

tY"~' tY"-J 



A' B 

^^n\2, (3.13) 



r 



12 



where we defined n\r, as 



12 



ni2 = ^. (3.14) 

ri2l 
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Hence, we obtain 



dPl 17117712 



dr- .?/"■- P-15' 



At the Newtonian order, we have Q\ = from Eqs. ( |B1D and ( [B^ ) , so we have the 



Newtonian momentum velocity relation, P| = mAV^. Hence we obtain the Newtonian 
equation of motion for relativistic compact binaries. 

IV. THE FIRST POST-NEWTONIAN EQUATION OF MOTION 

In this section, we shall derive the IPN equation of motion. At this order T^^ is trans- 
formed as 

T^^ = (rA)'Tr + 2e\\TA + 0{e'), (4.1) 

where Ta is a transformation coefficient whose explicit expression is irrelevant for the purpose 
in this paper. The second term of the above equation does not contribute to the near zone 
dipole moment when we require that J\ vanishes, so the condition D\ = defines the same 
center of the mass as in the Newtonian case. However, we must take into account this 
transformation when we consider the spin-orbit coupling force as seen in the next section. 

As for Q\, this term appears from the 0.5PN order and depends on the size of the body 
zone. However up to the 0.5PN order, Q\ does not affect the equation of motion. We will 
show this explicitly in the following two subsections. 

A. The boundary dependent terms 

Before deriving the IPN equation of motion, we present some comments on the mass 
and momentum defined in this paper, which are relevant to the cancellation between terms 
depending on the size of the body zone; our definition depends on the size of the body zone 
since the pseudotensor has non-compact support. The boundary of the body zone is given 



12 



by hand so that the mass must depend on the size of the body zone. Thus, we may spht 
the mass as 

mA = rhA + ruA, (4.2) 

where rfiA does not depend on the boundary size but tua does. In order to perform this 
sphtting, we use the fact that the integral of A]^ over the whole space W does not depend on 
the size of the body zone, say eRi. In fact, from the relevant expression of the pseudotensor, 
the equations ( P^ ) in the appendix B, we obtain at the lowest order 



Jw/Bi Jw/Bi 



N 
Bi 

7 roo 






+ (terms independent of Ri) 

7 TTl^ 

= — e^ — -^ + (terms independent of Ri), 
2 eRi 

where W/Bi means the spatial domain obtained by subtracting Bi from W. Since the total 
energy must not depend on the size of the body zone, we find 

^A = e'l^^+0{e'). (4.3) 

In the similar manner, from Eq. ( [B5| ) we split P^ into P4 = P^ + P\ and obtain 

PA = e2-!^ + 0(e2). (4.4) 

^ 3 eRA 

Furthermore, we evaluate each integral in Eq. ( p^.25D as 

7ml 



e ' <p rf5,(6)A]v^{x^ - z\{r)} = -e'—^, (4.5) 



2 7 



e-^l rf^.(6)A]vH^^-4(r)} = -e^^. (4.6) 



IdB, ^''" '' ^ '' '^ eRi 

Hence, we obtain Q\ = O(e^) and 



2 i 

2^'^A^A , nit Ji\ 



Q'a = ^'^ + 0{^\ (4.7) 
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where we split Q\ like mA- Thus we have 

P\ = mAv'j, + 0{^). (4.8) 

This is the momentum velocity relation independent of the size of the boundary. Note that 
this relation is valid up to the 0.5PN order, since Q\ and all the terms which depend on the 
size of the body zone are evaluated with the IPN pseudotensor. There are other terms which 
depend on the size of the body zone. These come from the integrals which are formally of 
the 2PN and 2.5PN order (e.g. (8)Ay in Eq. ( [4.5|) ) and, at first sight, contribute to the 
IPN equation of motion. We will discuss this issue in the forthcoming paper where we will 
derive the 2.5PN equation of motion. 

Since the boundary-dependent term of the mass is of order e, when we pay attention 
only to the Newtonian order, we can rewrite the Newtonian equation of motion as 

,-„,M = _!!!l!!!i„;,. (4.9) 



dr r 



B. The first post-Newtonian equation of motion 

In this subsection we derive the IPN equation of motion paying attention to the cancel- 
lation between terms dependent on the size of the body zone. 

First we consider the temporal component of the four momentum at the IPN order. 



From Eqs. O, (g|), (|3J) and (gj), we obtain 



CLl' J oB\ J oB\ 



217111712 



(4(fi2-t;i)-3(fi2-^2)), (4.10) 



' 12 

where we used the following integral formulae 



IGvr 



JdB^ 3ri2 3rf2 



;rr,i Tri,k ^mim2v{8'^r\2 , Amim2vi5l^rl2 



rf^m(4)/i (4)ft = 5-3 + 5-3 • (4.12) 
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Using the Newtonian equation of motion, this is rewritten as 



dP[_ 
dr 



ar^2 ri2 



(4.13) 



Hence, we obtain 



PI = mi 



a/1 n Sm^ 



+ 0(e' 



(4.14) 



-2 - ri2 

where we renormahzed the integration constant into rrii and replaced rfii by ?7ii at the e^ 
order because the difference is of order e higher. We spht P^ into the boundary independent 
part PJ and the boundary dependent part PJ. Up to the 0.5PN order P^ = niA- 

Now we turn to the equation of motion. The relevant expressions for the pseudotensor 
and the metric components are given in the appendix B. The surface integrals appearing 
here are as follows p2| ; 

JDT JDT 

1 2 „ 



(9B] 



dBi 



dSk{i)y—gt\ 
dSk(e)[-gt\ 



2 '^121 

ri2 
rhiTfi2 i 

2 ^^12 

ri2 



(4.15) 



■-{v2-nu) +-V2 



8 ^ 8mi m2\ 

3 ri2 ri2^ 



77111712,^ ^ . ,• 77111712 

{v2-ni2)vi + —^^ 

'19 J 



^12 



12 



(J 
(-(^^1 -^12) -3(^2- ^12))t'2 



2llmjm2 



n 



12 5 



t/S", 



asi 



fc(6) 



-gt 



LL 



3ei?irf2 

—^^12 I 0(^1 -^2) 



^12 
4:7hi7h2 



77117712 



[V2 ■ 7li2)Vi 



3r2 



[7ii2 ■ t;ijt;2- 



(4.16) 



(4.17) 
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We note that the last term in Eq. ( [4.16| ) precisely cancels out the time derivative of P\, Eq. 
( [4.4| ) (or equivalently, m^f^ + Q\) when we use the Newtonian equation of motion ( ^4.9| ). 
Inserting these results and Eq. ( [4. 141) into Eq. (|2.27D and using Eq. ( [4.8|) , we obtain 

7hi7h2 



dv\ 

7711-— 

dr 



2 ^12 



+e 



ri2 

27hi7h2 



r 



12 



^12 1 - 



/ 2 „2 3,^ ^,9 ,.^ ^x 5mi 4m2\ 

-i;i' - 2^72 + 7: t?2 ■ ni2f + 41/1 ■ i;2 + - + -) 

^ 2 ri2 ri2 ^ 



+ F*(4(t;i-ni2)- 3(172 -^12)) 



(4.18) 



where we defined the relative velocity as 
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V = Vi- V2- 



(4.19) 



Thus, it has been shown that the terms dependent on the size of the body zone do not 
affect the 1 PN equation of motion, since parts proportional to R^^ cancel out totally. 

The equation (^4.18|) takes the same form as those in previous works |19,11|. Hence, the 
applicability of the IPN equation of motion has been extended to compact objects. 



V. THE EQUATION OF MOTION WITH EFFECTS OF THE QUADRUPOLE 

MOMENT AND THE SPIN 

Now we consider the forces associated with the spins and quadrupole moments of the 
component stars in the equation of motion which appear at the 2PN order in our ordering. 

First, we must derive the metric components which depend on the spin and the 
quadrupole moment. The calculations are straightforward except for the spin part of the 
field (8)^^^- From the equation (|4.1|) and the requirement that J\ vanishes, we can calculate 
this term as 



h^ 



4e^ ^ / d'aA- 



A]J{t - t\rA - e^dA\,aA] e) 



+ 4 



\rA-e^aA\ 



a y 



4e' E 

A=l,2 \JA JBa 



— / rf^aAe^A]v^(r, a^; e) + ^ ^^ 



r" 



Ba 



ic/B \x — y 

d^aAe^a^A^^ir, Sa; e) 



+ (terms irrelevant to the spin and higher order terms than 2PN) 
va r , 



4^^ E 



A=l,2 L 



:-A + e"^{^d\ + ^Mfv\) 



+ (terms irrelevant to the spin and higher order terms than 2PN) , 



(5.1) 



where in the last equality, we used the equation (2.19) in [^]. We defined d\ and the spin 
tensor M^^ of the body A as 



d\=\ime^ d^aATlahAV, 

c^O Jba 

Mf=\im2e' [ d'aAa^AJ,\ 

•^^O Jba 



(5.2) 
(5.3) 



Recalling that AY is a component of a tensor in the Fermi normal coordinate, we define the 
center of mass by requiring that d\, instead of D\, vanishes. 
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It is important to note that the above choice of the center of the mass leads us to the 
following velocity-momentum relation; 

Pl = Plv\ + Q\-e' '^'^^^'''\ (5.4) 

^12 



Here we used the fact that the temporal derivative of the spin M^ is of order e^ [24 
The required metric components are expressed as 

hi h 

A=l,2 'A A=l,2 'A 



A=l,2 ''^A 



(6)/.^^=4E 4mM 



A=l,2 ^A 

where J^ denotes the tracefree tensor for the reduced quadrupole moment of the body A 
defined as 

Pi = lirne^ / d^aA^Yia^aj^ - — -la^P). (5.5) 

The spin vector S*^ is related with M^ as 

S\ = le,,,Mi\ (5.6) 

where eijk means the spatially unit antisymmetric symbol. 

The pseudotensor needed in this section is expressed as the equations (|B8|) , ( p9|) and 
( P10|) in the appendix B. 



A. Quadrupole-orbital couplings 

First we calculate the contribution from quadrupole moments. The surface integrals that 
we need to evaluate are only 

3 

12 

From Eq. ( p.27|) , we obtain immediately the result: 



/ dS.^J-gtfr) = -7^(mii'2' + m2/f 1(2^^2 - 5n\,n1,n[X (5.7) 

JdBi ^ ' ^^\2 
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mi^ = e'^(mjf + m2/f ) {^^''n^, - Snl^nfX^)- (5-8) 

Here we omitted the Newtonian and the IPN force and the spin-orbit couphng force. Note 
that we use ffiA instead of PJ in the left hand side of the equation ( p.27D here since the 
difference between P\ and friA does not affect the quadrupole-orbital couphng. This is 
because from the equations ( pSf ) and (|B9D the equation of evolution for time component P\, 
Eq. ( p.26| ), does not contain the quadrupole moment at the 2PN order, hence P\ does not 
contain the quadrupole moment. 

This equation, (|5.8| ), has the same form as the Newtonian equation of motion for the 
so-called tidal force. Therefore, it has been shown that the above equation is valid even for 
the compact objects. 

B. spin-orbital couplings 

Next we turn our attention to the spin-orbital coupling in the equation of motion. The 
relevant surface integrals are as follows: 

-3 2^1,,,-. 2mi 



£^ ^^.(8) (-.tiy = e^l^Mf - |^(M^'A'^' + M-A' 



2m 



ISrfs 



(MfA'^' + 3MfA'^), (5.9) 






m2 



3-Ri ' 12 

6 



r 



12 





+ v\[--M^A^^ - 2Mf A'*)] , (5.10) 



where we defined 



A^^' = 5*^' - Sniania- (5.11) 

We also need to evaluate the surface integrals in Q\ which are 

^^^ d5,(8)(-(7tiy {x^ - zi(r)} = ^Mf <2. (5.12) 



Then we obtain 

Q\ = e^^Mf <2. (5.13) 



,4 ^"^2 j^jik^k 
■J' 12 



Note that the surface integral of (8)(— 5'^ll) does not depend on the spins of the bodies. 

Inserting the above results into Eq. ( |2.27| ) and recalling the velocity-momentum relation 
(E^, we obtain 



dv\ ,V^ 

mi-— = -e ^T- 

dr rf2 



(2miM| + maMf ) A'^ + 2(miM^^ + msMf ) A'^] , (5.14) 



where we omitted irrelevant terms to the spin-orbit coupling force. Note that the difference 
between P^ and rfiA has no influence on the spin-orbital coupling because of the same reason 
as in the quadrupole-orbital coupling. 

Here, it is worthwhile to mention that the terms proportional to -Rf ^ cancel out totally 
in the equation of motion as we expect. This equation is exactly the same as the equation 



in pl|22 



VI. CONCLUSION 

Based on the point particle limit, we derived the equation of motion not only at the 
first post-Newtonian (IPN) order but also with effects of the quadrupole moment and the 
spin of compact objects at the 2PN order. It is worthwhile to mention why these effects 
are of the 2PN order, though they are usually of the Newtonian and the IPN order. In the 
point particle limit, the ratio of the size of compact objects to the separation of the binary 
is taken as an expansion parameter, which is equal to e^. Hence, the quadrupole moment 
and the spin become of higher orders in the ratio of the radius to the separation. Also, we 
assume that the velocity of the spinning motion of the body scales as e. Therefore, their 
effects appear at the 2PN order. 

As for the spin-orbit coupling, we must carefully define the center of the mass of the 
object to have the same expression as the previous works ||T^^,^. It is well-known that 
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the form of the spin-orbit couphng force depends on the definition of the center of the mass 
of the object p8 |. 



We emphasize that the definitions of the mass, spin, quadrupole moment of objects 
include the effect of the strong self gravity explicitly (for the case of the weak internal gravity, 
e.g. |2^J5^). Also we note that we did not need any regularization schemes. As a result, we 



have shown that the post-Newtonian equation of motion can be applied to compact binaries 
with strong internal gravity when we use the suitably defined mass, spin and quadrupole 
moment up to the order studied in this paper. We are pushing the present scheme ahead to 
complete the 2nd and 2.5 post-Newtonian equation of motion for the strong self-gravitating 
objects which will be investigated in the near future. It is also very interesting to see if the 
present scheme works for the 3rd post-Newtonian order. We hope to tackle this problem 
also in the near future. 
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APPENDIX A: COMMENTS ON THE DERIVATION OF THE METRIC 

The mass and four momentum are defined as the volume integral of A^'^ in the body zone 
so that these depend on the size of the body zone. On the other hand the integrated form 
of the reduced Einstein equation ( p.l9[ ) has no information of the size of the body zone, so 



we expect that the metric components do not depend on it. We here show that the field h'^ 
is really boundary-independent up to the 0.5PN order. 
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From Eq. ( |2.19| ), we have up to the IPN order 

yltt,2 '^A Jc/B \x-y\ 

We can evaluate exphcitly the boundary-dependent term arising from the non-compact 
source, i.e. the last term in the above equation. 
From the Eqs. (p^ ) and (|3.3|) , 

^3^ (6)Ay(^.^') = _JL V / Sy , ^ , ""i ^, + (irrelevant terms). (A2) 

c/B \x-y\ 87r^^_2-^c/i3 |?/^|^|a;-y| 

Then it is easily verified that 

d^y . ^ Ml ^ ^= / c?^l/i I ^ Ml ^ ^~r + / ^^Vi I ^ Ml ^ —- + (higer order terms) 

At: f I 1 \ 

= — — + (higer order terms). (A3) 

ri Veiti 2ri/ 

Inserting the above result into Eq. (^1|) we obtain 



,...4.«i:;^(i'i-^||)+0(a. (A4) 

Let us split P\ into P\ and P\ as P4 = PJ + PJ where P\ depends on the size of 
the boundary while the P\ does not. Then from the equation ( [4.14| ) we have relations like 
P\ = rriA, P\ = rfiA and P\ = tjia up to the 0.5PN order. In section IV we obtained the 
following relation 

rfiA 



,7m\ 



2eRA 
Thus it is shown that up to the 0.5PN order the boundary-dependent term in P^, i.e. P^ 

cancels out the term arising from the non-compact source. We can write the metric up to 

the 0.5PN order as 

/i- = 4e^y:^ + 0(e6). (A5) 
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APPENDIX B: THE METRIC AND THE PSEUDOTENSOR EXPANDED IN e 

For convenience, we collect some formulae here. First, we show the components of the 
Landau-Lifshitz pseudotensor expanded in e up to an order relevant in this paper. The 
lowest order in the pseudotensor is e^ since the lowest order of the near zone field h'^'^ is e^ as 
shown later. The components at the fifth order are identically zero, while the components at 
the seventh order vanish because {5)h'^'^ and (7)/?-^^ depend only on time r ||2^. The required 



orders of the gravitational field and the Landau-Lifshitz pseudotensor to obtain a certain 
order equation of motion are summarized in the Tables || and |T|. 

( ) and [ ] mean the symmetrization and the antisymmetrization respectively. 

O(e^) 



(4)[-167rc/t^\] 



0, 

0, 
1 



l,TT,i l,TT,j 

[4)11 (4)ll 



-S''\4)h'^'^' {4)h'^'^,k 



0{e') 



(6) 
(6) 
(6) 



■IGngf 



7 



LL\ 



l,TT,k -ITT 

(4)11 {4)11 ,fc 



[-167r<V] 



2(4\h'^'^ h(A\h7 '*' -(4\h'^ k(A\h7'^' 



'.{4)11' ,fc(4) 



4 



{4)11' ,k{4)l 



-IGngq^] = 4{4)h-^''''\4)h\k''' - 2{4)h7^^\4)h''\k 



+ 2(4)^ " (6) 



1 - - 1 - - 

TT,(i U\tt\,J) I ^.TTXi U\k J) J^TT UTT,l UTT,1 

''-n.' + ^(4)'^ {4)1^' k\ - -{4)11 {4)11' (4\ri 



(4)/ 



:ti 



kl 



'TkA 



— -{4)h7' ,fc(4)^^ ,/ + {4)h'^'^ ,k{4)h ,/ + (4)/i^ ' {4)h'^lk 



l,TT h''"''^'^ h 



:{4)n ,fc(6) 



\TT,k ul 



1 



h i.k + -(4)h'^'^(4)h'' 



{4)11' ' (4)"' /,fe -^ 7(4)"' (4)"- ,fc(4) 



^h 



TT.k 



0{e') 



{r)[-16ngtll] 
{r)[-16ngtl\] 
{7)[-167igfl^] 



0, 
1 
2 



(4)/l (5)/!' kl'-" + {4)h (7)/!' 
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(Bl) 
(B2) 
(B3) 



(B4) 
(B5) 



(B6) 



1 ..„• 



-5'^ {(4)^^^' (5)^ k,l + (4)^^^' (7)^^^,fc} 



Ofe^ 



(8) 



—IQTTgt 



LL\ 



7 

'a' 



t,TT,k ITT 

r(4)«- (6)/^ ,k 



(4)/i^ ,fc(4)^^ ,« + (4)^^ ' {4)h'^{k,l) 



\kl 



+ -(i)h'^^' (4)/l /,fc + {4)h'''^ ,k{i)h ,1 



(8)[ 



■167r5't 



1 



llJ 



'2{4)h'^ k,l{4)h ''*' ' + {4)h'^\k{4)h ,/ + T(4)/i^ ,/(4)/i 'fc* 

1 



4 



(4)/i^^'*(4)^ /,r 



+ 2(6)/l^"='^l(4)/l 



-(p.^hT' k(A\h 



4 



(6)'i ,fc(4)' 



+ 2(6)/i^'^,fc(4)^^' '*' — T(6)^^'^'\4)^^ ,fc 

3 - 

+ 2(4)/l^^(4)/l^^^/(4)/l^'*'' + -(4)/l^^(4)/l^'^'*(4)/l^ i 
~ 7(4)^'^ {4)h'^'' ,k{4)h'^'^'' + -(4)^^* (4)^^^,/ (4)^^^' , 

4 o 



(8)[ 



■16'Kgf-' 



LL\ 



-2(4)/l^-'<\6)/^^'',Z + -(4)/i^^'^6)/^'-^| 



' ,.i) 



(B7) 



[B8) 



(B9) 



+ 2(6)/i^'''\4)/i^'\/ - 2(6)/l^(*'l'l(4)/l^'^^ + 2(4)/l^''<\6)/^^'^\/ - 2(6)/l^z'^\4)/il""'^'^ 



;(4)'i (8) 



+ 5'^ 



— -(4)^^ ,fc(6)^'^ ,/ + 2(4)/l'^ ' {6)h'^[k,l] 



\TT,k ul 



+ {4)h'^^ ,l(6)h ,fc — T(4)^^ ' (6)"- «,fc 



/l' 



l,TT,K -ITT 

:{4)ii. {8)/^ ,fc 



+ (Terms independent of the spin and quadrupole moment) . 



(BIO) 



Next, we show here the field h^'^ expanded in e. The derivations can be seen in |^^, ||17|| , 
[151 and partly in the section V and appendix A. We do not show the field of order e^ and 



e^, since these depend only on time ||2^ and hence have no effect on the equation of motion 
up to the order relevant to this paper fS^ . 
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/,- = e\,)h^^ + e\e)h^^ + e\^)¥^ + e\^)}f^ + 0{e^), (Bll) 

(4)/.^- = 4 E ^, (B12) 

A=l,2 ^A 
(6)ft =-2 2^ U^^A-^a) -t^Ai + S 2 ^12-(«l-«2) 



A=l,2 



^A ^12 



r/i^ 17111712 mim2 



+ 7 V ^ + 14^^^i^ - M^^^i^ V -, (B13) 



A=l,2 'A A=l,2 'A 



+ (Terms independent of the spin and quadrupole moment) . (B14) 



/i- = e\^)h^' + e\e)h" + O(e^), (B15) 

(4)/^" = 4 ;^ ^^^, (B16) 



A=l,2 



Ta 



r 



k 

A ■\fki 



A=l,2 'A 

+ (Terms independent of the spin and quadrupole moment). (^17) 



h'^ = e\^{h'^ + e^(5)P + €"^(6)^ + O(e'), (B18) 



(4)/^'^' = 4 E 

A=l,2 ^^ 



+ E 



^2 



rriA i j ^mim2 ^ j 

A=l,2 ^A ^12'^ 



-8 

(6)h'^ = 4 



(^'fc^^^ - ^S'^M 


mim2,^ 

g2 (^12- 


- niY\ni2 + n2y\ 


k 
\ ' 'A ]\/Tk{t j) 

4=1,2 ^A 







(B19) 



+ (Termsindependent of the spin and quadrupole moment). (B20) 
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Near zone 



Near zone 



Body zone 



[T,a- 



rptt 



T*; ~ e- 



n - ^ 



~ e 



A ~e 



n - e- 



rri^ —8 



TABLE I. Scalings of the stress energy tensor T'^ of matters on the initial hypersurface in 
various coordinates. 





Newtonian order 


1 PN order 


2 PN order 


h7^ 


e^ 


e6 


e« 


/i" 


— 


e^ 


e6 


h^i 


— 


e^ 


e6 



TABLE IL Required order of the gravitational field to obtain an equation of motion (EOM) up 
to a certain order. For example, to obtain the Newtonian EOM, we have to calculate h7'^ up to e^. 





Newtonian order 


1 PN order 


2 PN order 


J.TT 




e^ 


e« 


fTi 


— 


e^ 


e« 




e" 


e^ 


e« 



TABLE in. Required order of the Landau-Lifshitz pseudotensor to obtain an equation of motion 
(EOM) up to a certain order. For example, to obtain the Newtonian EOM, we have to calculate 
*LL up to e^. 
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